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Abstract. We show that superradiant emission (SR) in semiconductor laser diode structures 
is governed by the master equation of the form of Ginzburg-Landau equation or Gross- 
Pitaevskii equation for a phase transition to coherent matter state. We conjectured 
condensation of one-dimensional longitudinal lower polaritons with effective mass of 
10~ n the mass of free electron. Two different regimes of SR emission and polariton 
condensation are predicted, one of which exhibits some characteristic features of the Dicke 
superradiance in atomic gases. The predictions of analytic model are confirmed by the 
results of numerical simulations based on Maxwell-Bloch equations, utilizing InGaN/GaN 
heterostructure quantum well as a model system. 
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1. Introduction 

The search for spontaneous build-up of macroscopic coherences in semiconductor 
microcavities has resulted in intensive research towards a polariton laser. Interestingly, 
group III - nitride semiconductors are capable of reaching macroscopic coherences at 
room temperature conditions. Thus due to large reduced exciton mass, the InGaN/GaN 
quantum well (QW) exciton binding energy is higher as compared to GaAs/AlGaAs QWs 
or other conventional III-V counterparts. Furthermore, due to a stronger exciton-photon Rabi 
coupling, InGaN/GaN microcavities offer possibility of polariton lasing at room temperature 
conditions (Q. High critical temperature of two dimensional (2D) BEC transition is 
conditioned by low in-plane effective mass of polaritons, so as the critical density is reached 
before destruction of excitons. Such 2D BEC is a transient dynamic state, which can 
nevertheless be analyzed using a stationary Schrodinger equation J2]|. 

The superradiance (SR) in semiconductor edge emitters can be considered as another 
example of spontaneous macroscopic coherences in solids 0. The cooperative radiative 
recombination in an ensemble of quantum oscillators (e.g. atoms or molecules) has been 
predicted before the invention of lasers 0. Since that, it has been extensively studied both 
theoretically and experimentally J5l|6]|. The characteristic features of the SR emission are the 
temporal and spatial coherence, highly anisotropic emission pattern, quadratic dependence of 
pulse intensity on the number of excited atoms I oc n 2 , afterpulse ringing attributed to Rabi- 
type oscillations. The SR pulse duration decreases with the number of emitters r sr oc 1/n, 
while the pulse energy is proportional to the ensemble population (Ir sr oc n), in agreement 
with the energy conservation considerations. On the other side, the spontaneous nature of the 
transition to the transient macroscopically coherent state is responsible for large fluctuations 
in the shape, duration and amplitude of SR pulses due to quantum-mechanical uncertainties 

OTH. 

In semiconductors, a hypothesis has been drawn that the SR is assisted by formation of a 
transient BCS-like state of electron-hole (e-h) pairs mediated by photons [fTOl [TTTl . According 
to that, the e-h system undergoes a second order non-equilibrium phase transition when the 
coherent e-h BSC-like state is building up during SR pulse emission. 

In this paper, starting from semiclassical travelling wave Maxwell-Bloch equations, we 
show that in semiconductors, superradiance regime is governed by the Ginzburg-Landau 
equation (GLE) or Gross-Pitaevskii master equation for a phase transition to coherent matter 
state. The effective mass of condensing one-dimensional (ID) longitudinal polaritons is 
very low, 10~ n m , giving a good reasoning for the thermal de Broglie wavelength criteria 
for macroscopic quantum degeneracy being not impacting our GLE for SR. This allows 
us to define the critical density in the system, the order parameter and the coherence 
length without accounting for the system temperature in the first approximation. These 
arguments thus further support the hypothesis of BCS-like condensation during SR emission 
in semiconductors. Our GLE for SR in semiconductors has two types of solutions, indicating 
that two different regimes of SR emission and polariton condensation are possible. One of 
them (type-I SR) reveals features that are similar to the Dicke superradiance in atomic gases 
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with long natural decoherence time. In this regime, the macroscopic coherence is established 
over the entire sample. In the type-II SR, only partial coherence can be achieved in the 
sample so as the domain of condensate fraction is smaller than the size of the sample. SR 
pulse parameters predicted by analytic model are confirmed in numerical simulations based on 
the travelling wave Maxwell-Bloch equations. Inspired by remarkable features of the group 
III - nitride microcavities in reaching room temperature BEC macroscopic coherences, we 
conducted our numerical studies of the SR in an edge-emitting ridge- waveguide cavity with 
InGaN/GaN QWs. 

A serious study of cooperative coherent effects in atomic (molecular) vapor has been 
started in 1970s, 20 years after the predictions of Dicke [9]. These studies revealed that SR 
emission, generation in masers, gyrotrons are of the same nature. The SR in small-size and 
extended samples has been studied. The type-II SR regime in semiconductor edge emitting 
devices predicted here might be considered by someone to be the Dicke SR regime in atomic 
gasses in case of an extended, optically thick medium [fT2l . However the differences are (i) 
in the relationship between the decoherence time and SR pulse width and (ii) the presence of 
a waveguide in semiconductor edge-emitting cavity. As a result, the Dicke SR in extended 
atomic gasses reveals the pulse intensity oc n 2 flU, while our predictions for the pulse intensity 
in type-II SR regime in semiconductor cavities is oc n 3 ^ 2 , the result which we confirm by 
numeric simulations. 

Because of low effective mass of longitudinal polaritons, we concluded that thermal 
potential would provide a negligible correction to the critical density for SR emission and 
polariton condensation. In 70s, there was quite large number of sophisticated calculations 
for SR phase transition in atomic or molecular gasses that specifically look at the critical 
temperature of the transition [fT3l . Following Dicke, they treat an ensemble of atoms (or 
molecules) as a correlated spin states lfl4ll . Unfortunately, most of these theoretical works are 
limited to obtaining a theoretical expression for the critical temperature T c , without attempting 
to make an estimate. Following along the lines of these studies, (e.g., examining Eq.(32) from 
Ref. |[T5l ). we conclude that the thermal potential can be accounted for in our expression for 
the critical density (fT3l) as 

r go K r -n t ) ih, 

v g (ai+2/L*) l 2k B T c ; ' 

The correction term from the thermal potential in the r.h.s of this expression oc 2 exp(^^-) is 
negligible, about 10~ 51 . 

The thermal potential k B T does not enters explicitly in threshold conditions for laser 
or maser equations, and this well agree with experimental observations. In Ref. lfT6ll . a 
possibility to derive GLE for a laser has been evoked. This conjecture is inline with our 
findings for the SR GLE (fT5l) that does not include thermal potential k B T. 

The paper is organized as follows. In Sec. [2l we introduce traveling wave Maxwell-Bloch 
equations used as a starting point of our development and in the numeric simulations. In Sec. 
|3j we obtain our GLE for SR and longitudinal polariton condensation in semiconductors. We 
analyze its two possible solutions. In Sec. |4l we provide a detailed discussion to our findings. 



Longitudinal polariton condensation and superradiant emission 



4 



2. Model system 

The SR in semiconductors is usually achieved in an edge emitting laser diode cavity 
configuration ifTTll . Such cavity has a waveguide structure directing the emission of 
spontaneous photons and lasing modes. Herewith we also consider the built up of SR in a 
semiconductor cavity pumped high above the transparency conditions of the semiconductor 
gain medium (bulk or quantum wells separate confinement hetero structure). Thus if during 
the build-up of SR pulse, a phase transition to macroscopic quantum state takes place, the 
condensation is essentially one-dimensional. 

In practice, to prevent the onset of lasing emission, a tandem cavity configuration with 
an additional absorber section is used ifTOl . The condensation takes place in the gain section, 
while the absorber section is just an auxiliary arrangement for experimental observation. It 
does not impact the physics of SR emission in strongly pumped semiconductor gain medium. 
Therefore, we limit development of our analytic theory to the processes in the amplifying 
semiconductor medium. 

As a starting point, we use a semi-classical Maxwell-Bloch equations for the traveling 
wave amplitudes, carrier densities and macroscopic coherences (polarization) in the ensemble 
of electrons and holes. (Similar equations can be found in studies of SR in atomic ensemble, 
in the case of extended source [12].) The evolution of the amplitudes of the forward (A + ) 
and backward (A J) travelling waves is defined by the cavity loss and macroscopic medium 
polarization 



where the field amplitudes are normalized in the secondary quantization convention 



The same normalization used for coherences (the ensemble average for the off-diagonal 
density matrix elements) induced between the electrons and holes and associated with the 
forward and backward waves reads 



Here the z axis is directed along the cavity waveguide. As shown below, the macroscopic 
variables P± provide a measure for the order parameter of the system. The dynamics of 
carrier populations and coherences is described as follows: 




(2) 



a±\N ± )=A±\N ± -l): 




(3) 




(4) 




(6) 



(5) 
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Table 1. Model parameters for wide-bandgap semiconductors, using data from Ref. lfl8l [T9l 
EO). 



Parameter Name Value 





QW width 


3 nm 


n q w 


number of QWs 


2 


Tig 


group refractive index c/v g 


3.5 


w 


waveguide width 


2 /im 


r 


optical confinement factor 


0.02 


n t 


carrier density at transparency 


1.6-10 19 cm" 3 


T n 


carrier lifetime (gain section) 


0.9 ns 




carrier decoherence time 


100 fs 


g Q =dg/dn 


differential material gain 


1.2-10- 6 cm 3 /s 



where the Langevin force term due to polarization (thermal) noise A± triggers the spontaneous 
build-up of superradiant pulse. This term is important for numerical solutions, while 
semiclassical analytical theory we develop here does not require it. 

In Eqs.©-©, v g = c/n g is the group velocity of modes in the cavity waveguide, g Q 
is the differential gain coefficient, r n and T 2 are the carrier relaxation time and decoherence 
time, the coefficient T accounts for the partial overlap of the waveguide mode field and the 
gain medium (e.g. quantum wells), n t is the carrier density at transparency. The coefficient 
cti accounts for the optical loss in semiconductor medium, J(z, t) is the pump current density 
and d is the thickness of the active region (e.g. d = Aqv/Nqw in a hetero structure with Nqw 
QWs of thickness dqw)- We neglect the carrier diffusion since it is not significant at the time 
scale set by the width of the SR pulse (~ 1 ps). 

For spontaneous emission, the coherence (polarization of the medium) decays faster 
then the population inversion. In that case, the medium polarization follows adiabatically 
the optical field (dP±/dt, dP±/dz = 0). The equation for polarization can be adiabatically 
excluded so as the model transforms into the standard travelling wave model for ultrafast 
semiconductor lasers. However in the case of SR, this adiabatic approximation is not valid. 

Thus, we use the same set of approximations and parameters as used in the traveling 
wave rate equation model of a semiconductor laser but without adiabatic approximation for 
the macroscopic polarization. In particular, the details of the vertical composition of epitaxial 
layers do not enters the model directly. 

In what follows, the model calculations are performed assuming InGaN/GaN double 
quantum well separate confinement heterostructure with the usual arrangement of epitaxial 
layers for an edge emitting laser lfT9ll2T| . The main parameters of the model are summarized 
in TableEJ 
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Figure 1. (a) Energies of longitudinal ID polaritons. The inset shows the band anticrossing 
at the band slope of v g /2. (b) Same energy diagram in the moving at v g /2 coordinate system. 
The dot indicates the LP polariton trap. 



3. Analytic model 

3.1. Ginzburg- Landau master equation 

A hypothesis has been previously stated ifTOlfTTfl that the SR in a semiconductor is mediated by 
a formation of the transient BCS-like state of e-h pairs. In order to interpret the SR emission in 
terms of BCS-like transition one should transform Eqs.© and © to the form of the Ginzburg- 
Landau equation (GLE) for the second-order BCS phase transition. The difficulty is caused 
by the fact the GLE is essentially a steady-state equation for an isolated system, while the SR 
is a transient process. It appears that to overcome this difficulty we have to change variables to 
the internal coordinates of the SR pulse that is, to consider the process in a moving coordinate 
system. Prior to reporting the technical details of such coordinate transform, it is useful to 
provide a physical insight into the nature of condensing quasiparticles and the meaning of the 
coordinate transform into a moving coordinate system. 

Figured] (a) shows the dispersion of the reduced energy of e-h pair in a direct band gap 
semiconductor in function of the wavevector component along the cavity axis. (The e-h pair 
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fulfills the energy and momentum conservation under the interband transition with emission of 
a photon). The e-h pair energy dispersion is plotted alongside with the dispersion of photons 
in the guided mode of the cavity. The scale is set by the wave number k of photons at the 
interband transition energy and is so small that the dispersion of e-h pair energy oc k 2 is not 
seen in Fig. [Q Since we are interested in a small region around k , the effects of optical 
waveguide mode cut-off or any other features yielding a deviation from the linear dispersion 
are not considered here. 

The coupling of e-h pairs with photons results in formation of lower polariton (LP) 
and upper polariton (UP) dispersion branches, which are the eigen solutions of the Jaynes- 
Cummings model [|22l l23l . We thus focus here on features related to ID dispersion (and 
condensation) of such longitudinal polaritons in an edge emitting cavity, while two other 
degrees of freedom has minor impact on the results reported here. 

Jaynes-Cummings model predicts that LP and UP band anticrossing occurs at k . For a 
typical photon number in the SR pulse (see Secj4]), the Jaynes-Cummings model predicts a gap 
of 5 meV (see the inset in FigfTJ). The group velocity of UP and LP at the band anticrossing 
is dE/hdk = Vg/2, a half of the cavity mode group velocity, indicating that longitudinal ID 
polaritons in edge-emitting cavities are the half matter - half light composite bosons. In order 
to meet the general considerations admitting ID condensation of polaritons, we shall now 
identify the energy trap for such longitudinal polaritons. 

Figure [H(b) shows the UP and LP band anticrossing plotted in the coordinate system co- 
propagating with UP and LP at v g /2 velocity. The relativistic correction at such coordinate 
transform is negligible. Thus l/y/(l — v 2 /Ac 2 ) ~ 1.01 for parameters shown in Tabled] The 
UP (LP) energies are obtained using the transformation E' = E — \hkv g . 

In Fig. [11(b), the band edge of LP with negative effective mass represents an energy trap, 
thus admitting condensation of lower polaritons. The effective mass of LP polaritons \rriLp\ = 
\dE /h 2 dk 2 \ ~ 10 _11 m assumes possibility of room temperature condensation, during SR 
pulse emission (m is the free electron mass). The thermal de Broglie wavelength criteria 
for macroscopic quantum degeneracy predicts very low critical densities, indicating that 
thermodynamic considerations do not provide the main limiting factor for LP condensation 
and SR pulse emission. Therefore it would be reasonable to expect that in the first 
approximation, the temperature will not enter the GLE for SR emission in semiconductors. 

The condensate, as we have shown above, is expected to move at v g /2 velocity. 
Technically, it is much more simple to perform analysis in a coordinate system co-moving 
with the building SR pulse in the cavity, as compare to the system moving at v g /2. The 
coordinate transformation to the internal pulse coordinates has received a wide spread 
use in the analysis of cooperative effects and resonance amplification in masers [|24l l25l . 
For simplicity, we consider only the forward traveling pulse, introducing new coordinates 
(=t — z/vg and z (instead of t and z, respectively). As shown above, in the case of SR in 
semiconductors, we may ignore the relativistic corrections. 

The phase transition under the question is related to the change of the effective phase 
relaxation time T| ff , which is comparable to the lifetime of photons in the cavity (l/t> 9 aj=3ps). 
The decay of macroscopic polarization due to dephasing of e-h pairs cannot be neglected. 
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On the other hand, the carrier injection into the active region (e.g. into QWs) and carrier 
relaxation at the rate l/r n can be neglected at the time scale of the pulse width, yielding a 
system of equations that is similar to the one used in studies of SR in atomic ensembles 

(7) 
(8) 

T . y ^-(n-n t )A + . (9) 

J-2 V J 2 

Usually, such system of equations has been solved by introducing an instantaneous (partial) 
pulse area variable, which transforms it into an equation for a non-harmonic pendulum 
|[T2ll24l . The anaclitic solution has been possible so far only for a point source considered by 
Dicke. In the case of elongated (extended) source, it has been solved only numerically. 

The key to our analytic solution of this system was an observation that the optical field 
of the SR pulse follows the macroscopic polarization of the medium within a certain domain 
of size L*. This can be seen from the features spotted by the numerical solution in Fig. [2] 
(Compare Figs, (b) and (d) as well as (c) and (e); further discussion can be found in Sec. HI). 
Therefore the analytical solution is obtained by substituting 1/L* — > d/dz with L* being the 
effective size of the domain with large macroscopic polarization of the medium (condensate 
fraction), yielding 

A + ((, z) = M L— P+ te, z). (10) 
V T 2 Vg{ai+2/L*) 

This original anzats allows us to exclude the wave amplitude from Eq.© and convert our 
problem of SR in a cavity waveguide into a task analytically similar to the well studied case 
of SR emission from a point source: 



dP + = Tg (n-n* cr ) 
d( T 2 v g ( ai +2/L* 
dn P. dP+ 



P+, (ID 

(12) 



d( (n-n* r ) d( ' 

where we have introduced a new parameter 

, v g ( Qi +2/L*) 
n cr =n t +^ — , (13) 

r# 

which, as shown later, is the critical density for the condensation phase transition. Following 
recipes established for a point source, the last equation in (fT2l) allows one to recover the Bloch 
vector conservation in a cavity waveguide 

(n(C, z)-< r ) 2 +p2(C, z)=(no-rC) 2 (14) 

with n = n(—oo) being the initial carrier density in the system, before the emission of SR 
pulse. 
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However, taking the second derivative in the first equation (fT2l) and substituting dn/d( 
and dP + / d( from (fT2l) . we obtain the following master equation for SR emission and polariton 
condensation in an elongated semiconductor sample (cavity waveguide): 



d 2 P ++ T 2 g 2 L* 2 



(n -n* r f-2P 2 



P+=0. (15) 



dC 'T^ g (2+ ai L*: 

This equation has a canonical form of the Ginzburg-Landau equation (GLE) (or Gross- 
Pitaevskii equation), in which P + is the order parameter of the system and the internal 
coordinate (=t—z/v g plays the role of spatial variable. For the carrier densities above critical 
n* CT (see Eq.(fT2l)). its steady state solution defines the condensate fraction. The analytical 
similarities between Eq. (fT5l) and GLE, allow one to assume that for n>n* r , the evolution of 
the order parameter is defined by hyperbolic secant function: 

P + (C,^)=Posech(C/r p ) (16) 

Substituting in (fT5l) . we obtain the parameters P and r p : 

P = n - n* cr , (17) 

(18) 



T 2 v g ( ai +2/L* 



Tg (n -n* cr ) 

At the peak of SR pulse, the order parameter P thus reaches the highest value admitted by the 
integral of motion (fl4l) . 

Using normalization convention © and relationship (flOl) we find the evolution of photon 
number in the cavity z) and at the peak of SR pulse 

A 2 At-z/v g )=-^ \ . (19) 

+ 9 g r 2 cos h 2 (i^) 

N -a 2 (Q)- Ta - r ^ o(r2 ° ~ (20) 
N+ - A+{0) - 90 T 2 -T 2 v 2 (a t + 2/L*y (20) 

For convenience and in order to shorten our formulas, we use the following convention for the 
pulse width and coherence time ll26l 

t c = 2r p (21) 

The pulse width r c oc l/(n — n* cr ) and the peak power density hwv g N + oc (n — n* cr ) 2 exhibit 
the expected behavior known from the studies of SR in an atomic ensemble flU|7]|. 

We thus have obtained that SR in a semiconductor follows the GLE (IT3T ) for the 
phase transition to macroscopic coherent state and the medium polarization P± is the order 
parameter of the system. The proposed consideration of condensation during SR emission 
is different from all previous treatments. It does not involve the system temperature. As 
shown above, the effective mass of condensing quasiparticles (ID lower polaritons) is so 
small (~ 10 _11 m ) that the critical density is very low. This renders the approximation (fT5l) 
to be valid even though the thermodynamic temperature of the system is not accounted for. 
The thermal de Broglie wavelength criteria for macroscopic quantum degeneracy is satisfied 
at the carrier densities much below the one stated in Eq.(fT3l). 
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3.2. Type-I superradiance regime 



The features of the SR pulse (1201) are in perfect agrement with predictions of the Dicke theory 
if the coherence length L c =r c v g (as measured in the cavity) is longer than the sample size L 
(the cavity length ETlO . In this regime, refereed herewith as the type-I SR, the characteristic 
length of the condensate fraction L* is limited by the sample size L*=L<L C . From Eqs. (fT3l) . 



181) and ([201) . we obtain the critical carrier density Tier, pulse width parameter Tp and the 



(J) 



peak photon number N_ 

ng 1 ) >n > ni?=n t+ V ^l +2/L \ (22) 

1 9o 

m _ T 2 v g ( ai +2/L) 
Tp - rflbCno-ng?) ' ( } 

N (i)_ T 2 = r 2 (7o(no-^ ) ) 2 

The second critical density nh- 1 ' introduced here is defined in the next Section [3731 

The output pulse power can be obtained using the photon density in the cavity (l20l) and 
accounting for the reflection of the laser chip facets: 



p«(M)=Mi - R)v g n , 2 V t ! t0 _ g/c . (25) 

cosh ( — 

where i? is the reflection coefficient of the facet, Vt = w(IqwNq W /T is the mode cross-section 
area, w is the waveguide width, dqw and Nq W are the thickness and number of quantum wells 
and T being the optical confinement factor. The pulse envelope in Eq.d25l) accounts for free- 
space propagation of the pulse, and t is the delay to emission of the SR pulse. Note that 
within the semiclassical approach we use here, the time to emission of the pulse cannot be 
defined and requires quantum- mechanical consideration fl4][7]|. 



3.3. Type-II superradiance regime 

With increasing carrier density, the coherence length L c = 2t p v 3 reduces. At n > nir , 
it becomes shorter than the sample length L indicating that the size of the domain, which is 
occupied by the condensate fraction, is smaller than the sample size (L*=L C <L). We refer 
this regime as the type-II SR. As shown below, in distinguishing from type-I SR regime, the 
pulse peak power does not exhibit quadratic growth oc (no — n cr ) 2 . In principle, this is the 
only SR regime that can be observed in experiments with edge-emitting semiconductor laser 
cavities of several hundreds of microns length (see Sec. |4]). 

The critical density nir^ can be obtained from Eq. (|24)) . At this carrier density, the 
coherence length coincides with the sample size so as r c = L/v g . Substituting r c /2 in the 
l.h.s. of ([24]) and solving it with respect to no, we find the the critical carrier density for the 
type-II SR regime is 

^w^^(i+^)>e (26) 
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At n > ncr \ the effective size of the condensate fraction domain L* in (fT3l) and (fT8l 



is conditioned by the coherence length L*=L c =2v g Tp II " > and we obtain a quadratic equation 
with respect to Tp 11 ^ that reads 



r (ii)_ T 2 (^+l/7j /7) ) 

T5(o(^o-^i)-(^ 9 ai+l/^ //)x 
This equation has only one positive root, yielding 




l+aiV q T 2 



(28) 



(29) 



4r^ 22(no-nt) 2 v /r^ T 2 (n -n t ) 

In the type-II SR regime, the effective size of the condensate fraction domain is smaller 
than the sample size L c = 2v g Tp H ^ < L. At each moment of time, this domain is centered at 
the peak of the growing SR pulse. However, considering the build-up of SR pulse propagating 
through the sample, we shall split our sample into domains of size L c and sum up the 
contributions from all domains in the sample. In the laboratory coordinate system (t,z), 
we thus have to account for the fact that different sample domains do not contribute to the SR 
pulse simultaneously. 

However, these calculations are particularly straightforward in the internal coordinate 
system of the pulse (£, z) that we use here. We have to account only for incoherent 
superposition of wavelets emitted from different domains in the sample. Otherwise, if the 
domains are mutually coherent, the coherence length L c will not be smaller than the sample 
size L, leading to type-I SR regime. 

There are thus ~L/L C mutually incoherent domains. We introduce index j to enumerate 
these domains and corresponding radiated wavelets A + j. For the order parameter of each 
domain (polarization) and slowly-varying amplitude of the radiated wavelets, we use Eqs. 

_ in) 

p +J = n ° . "":„,< oo) 



cosKC/rf)) 
T 2 




(3D 



The random phases (f)j account for the mutual incoherences between domains of condensate 
fraction. The corresponding superposition of wavelets reads 

A + (C, z) = J2 A +tj e-^ {z - zUr) Q{z - z&), (32) 

3 

where the exponential factor accounts for material loss and 0(z) is the Heaviside step 
function, which accounts for the fact that a wavelet emitted at some z^> > z and propagating 
in the positive z-axis direction does not contribute to the field at z < z^ . 
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The superposition of incoherent wavelets (1311 ) defines the evolution of photon number in 
the cavity N + (£, z) = (|A|_(£, z)\ 2 ). In particular, the photon density at the cavity facet reads 

2\ 



2 



(33) 



jX e -«i(^-^ 0) ) (34) 



3 

T 2 1 - e- QsL 



,o(rr^l-e-^"cosh 2 (C/rW)- ^ 
Here the ensemble average is taken over the mutually incoherent domains of the sample. The 
cross-correlation terms between i-th and j-th domains vanish since the relative phase fa — <pj 
is random. Thus the factor (1 — e~ aiL ) /(l — e~ aiLc ) distinguishes the type-II SR from its type 
I counterpart (compare with Eq.(f24])). 

The peak photon number and the output pulse power in the type-II SR regime are 

N M = 1 ~ e ~ a ' L f36 ) 



p£?(m)=mi - r> 9 ^ (37) 



N 



(n) 



cosh 2 (^ 



(ii) 





where we have used that L c = 2v g r p . 

At large carrier densities n 3> riir \ the pulse width (coherence time) reduces, yielding 
the asymptotic behavior 

r™~i L_ „x , (38) 

iV| JJ) ~ ; ' (39) 

2g aiV g (r^ ') 3 

u -n t ) 3/2 (l - e- Qii ) (40) 

Thus in the type-II SR regime, because only partial coherence can be achieved over the sample 
(L c < L), the photon number does not follow the quadratic growth predicted by Dicke' model. 
At the same time, the growth of N+ is superlinear, as opposed to the linear growth in an 
incoherent ensemble. The pulse energy oc N+Tp is proportional to (n — n t ), in agreement 
with the energy conservation considerations. 

With increasing sample length L, the difference between the first n£) ((24)) and second 
nii 1 ^ (1261) critical densities reduces, indicating that the cavity tends to operate in the type-II SR 
regime provided special provisions are taken to prevent lasing. At very large sample length, 
both critical densities are clamped at n t + v g ai/Tg . It can be seen that some excess above 
the transparency carrier density is always required to compensate for material losses. 

If the optical losses are negligible (pa 1/L), the photon number and the output peak 
power in the type-I SR regime increase as ocL 2 provided the conditions of the type-I SR 
regime are fulfilled, in particular that L 2 <AT 2 v 2 /Tgo(n — n t ). 
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In the lossless limit of the type-II SR regime, the photon number and the output pulse 
peak power scale as cx L with the sample size. At same time if some losses are present in 
the cavity, the photon number and the peak power show a saturation, starting from the cavity 
length L ~ l/a*. 

4. Results and discussion 

The validity of analytical model predictions is confirmed by numerical simulations in an edge- 
emitting laser cavity incorporating InGaN/GaN double quantum well hetero structure (Table 
[Q). The numeric model used here is based on the traveling wave Maxwell-Bloch equations 
of Sec. [2] It reproduces spontaneous triggering of the SR emission by incorporation of the 
spontaneous polarization term A± in Eq.©. We also implement the boundary conditions at 
the left (z—0) and right (z—L) cavity facets in order to mimic the effect of accumulation 
of spontaneous photons in the cavity before the SR pulse emission and in order to analyze 
possible effects of multiple cavity roundtrips due to reflections of the cavity facets 

A+(0,t)=y/RA_(0,t), A_(L,t)=y/RA+(L,t). (41) 

We complete the numerical model by incorporating separate equations for a possible absorber 
section in the cavity. This enables us to compare the analytic model predictions for the SR 
pules and condensation of longitudinal polaritons (in the gain section) with the numerical 
simulations for the entire cavity representing the actual experimental arrangement. 

The model equations for the carrier population and polarization dynamics in the absorber 
section take into account (i) the Quantum Confined Stark Effect (QCSE) energy shift of the 
absorption edge and (ii) the reduced carrier lifetime in a negatively biased absorber: 




where a is the differential absorption to gain ratio, r a is the carrier lifetime in absorber section. 
The parameter riy accounts for the QCSE energy shift in the absorber section QWs caused 
by the external negative bias. In our numeric model simulations, we use a = 3.5 and r a = 5 
ps while the normalized absorber bias V a =(n t — n v )/n t is in the range form to -1 (highly 
absorbing state of absorber). 

Figure [2] displays the results of numerical simulations for the SR emission in the cavity 
of the overall length L=800/im incorporating a 160 /xm long absorber. The absorber section 
is located at the left facet of the structure (at z—G) and the gain section is situated at the right 
facet. The spatiotemporal dynamics of the carrier population (Figj2| (a)), order parameter 
(Figs|2tb),(c)) and the optical field (Figs. [2](d),(e)) as well as the intensities of the output SR 
pulses emitted at both facets of the structure (Figf3]) are quite complicated. 

Initially, a high reversed bias is applied to the absorber leading to depletion of carriers in 
the QWs (blue color in the scale of Figflta) ). It prevents the structure from lasing and enables 
accumulation of carriers in the gain section. The gain section QWs are pumped to a high level, 
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Figure 2. Numeric simulations: Spatiotemporal dynamics of the carrier population (a), 
macroscopic polarization in InGaN/GaN QWs associated with the forward FW (b) and 
backward BW (c) waves; field amplitudes of the FW (d) and BW (e) waves. The overall 
cavity length is 800 (z axis). The absorber of 160 /im long is situated at the beginning of 
the coordinates. The time interval of two cavity roundtrips is shown (t axis). All parameters 
are normalized on carrier density at transparency n t . The current amplitude is I = 281 mA 
(the current density J=22kA/cm 2 ). 
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Figure 3. Numerical simulations: Output SR pulses at the gain section facet (red curve) and 
absorber section facet (blue curve). The cavity length is 800 /im, the absorber length is 160 
/im and driving conditions are given in the caption of Fig. [2] 



so as in this particular case, the carrier population is about 15 times above the transparency n t 
(red color in FigJ^a). The entire process starts at t=0 from the spontaneous polarization noise. 
The rate of spontaneous polarization noise into the cavity mode ~^Ty/g T2A± corresponds to 
27 photons per cavity round trip. This cannot be seen in the scale of Figj2] The noise sources 
are uniformly distributed along the cavity but the macroscopic coherences are building up 
at the edges of the gain section, after about half of the cavity roundtrip time (T cav =18.7ps). 
Because of the asymmetric cavity, the shape and emission time of the output SR pulses from 
the left and right cavity facets are different (Figf3]) and the emitted SR pulses are not related 
by reflections at the cavity facets. 

The spatio-temporal dynamics of the order parameter P± shown in Figs. |2] (b) and 
(c) clearly indicates that during all process of SR, the coherence length remains shorter 
than the gain section length (L C <L S ). Interestingly, in Figj2] (e), the peak intensity of the 
pulse corresponds to photon density in the cavity of 1.4-10 19 cm~ 3 , while the density of the 
associated coherence excitations P± at the polarization peak maxima is of 9.9- 10 19 cm~ 3 , 
which is of the order of the initial carrier density n =2.1-10 20 cm -3 . Thus the associated 
photon density in the cavity is, at first sight, quite low. In fact the dynamic variables in Eq. 
© account for the optical confinement factor T due to the difference between the cavity 
mode size and the QW width. The effective peak photon density reduced to the QWs width 
is of 7-10 20 cm -3 , which is thus spectacularly large as compared with the initial density n 
of uniformly distributed carriers. In this example, illustrating a typical situation in the edge- 
emitting laser cavity, the system operates in the type-II SR regime, in which the condensate 
fraction is confined to a domain smaller than the size of the sample. 

The build-up of the macroscopic polarization is followed by Rabi oscillations at a 
frequency dependent on the pumping rate. They are clearly seen in the electric field 
amplitudes of the waves showing characteristic oscillations with the reversal of the sign 
(Figs|2d) and (e)). The carrier density drops abruptly at the same time when the polarization 
rises, indicating that almost all carriers contribute to the field (Figsffta)). When the travelling 
backward pulse hits the absorber section at z=160fim, the pulse intensity is sufficiently high 
to saturate absorber (the blue step in the region 0<2<160/im and t > 9ps in Figflta)). From 
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this time, there is no absorption and gain in the cavity. The backward travelling SR pulse 
reflected of the left cavity facet (at z=0) can now freely travel through the structure yielding 
emission of the secondary pulse at the right cavity facet at t ~ 19ps. Its amplitude can be 
appreciated from Figf2](d). (In the scale of Figf3l the intensity of this secondary output pulse 
cannot be seen.) In this particular realization of SR emission, the pulse emitted from the 
absorber section side is of smaller peak intensity. 

The peak power, FWHM pulsewidth of SR pulses and delay time to their emission 
vary within the range of 50-200 W , 0.4-5 ps and 9-36 ps , depending on the driving 
conditions (FiglU). Due to spontaneous polarization A± in Eq.©, the build-up of 
macroscopic polarization has a stochastic nature, which results in differences between 
individual realizations of SR pulses. In addition, large stochastic variations of the SR pulse 
FWHM width between adjacent data points in Fig@](b) are caused by the pulse ringing effect, 
when one of the satellite pulses, as in Fig. [3l reaches the half of the maximum of the main 
SR pulse. Therefore the lower boundary of data points in Figj4](b) should be regarded as an 
estimated width of the solitary SR pulse. For convince, we show the eye guides in Figs. H] (a) 
and (b) [black doted curves] . 

The time to emission (delay-time) is plotted in Fig@](c) as a function of the pump current 
density. At low pump rates, it takes around 1.5 cavity roundtrip until the emission of the SR 
pulse and condensation of lower polaritons. The SR pulses at two cavity facets are emitted 
almost simultaneously. At higher pump current, the SR pulses are emitted after approximately 
0.5 cavity roundtrip, following the scenario depicted in FiglU Note that higher pump levels 
shorten the time to emission of the SR pulses. 

Varying the overall cavity length and the relative length of absorber, we find that the 
intensity ratio and the width of the forward and backward SR pulses can be effectively altered. 
The common feature of all realizations is that the polarization and carrier excess above the 
transparency vanish after a few oscillations so as the the optical pulses travel freely without 
gain or absorption. This process is accompanied by building up of macroscopic coherent 
domain in the sample, which we attribute to condensation of ID longitudinal lower polaritons. 

Let us compare the results of numerical simulations with predictions of simplified 
analytical model based on solution of the GLE (fT5l) . In conditions of Figf2l the sample 
size is defined by the the length of the gain section (640 /im). The critical carrier density 
of n[r=2A8n t is achieved at the pump current density of 4.17 kA/cm 2 (no=T n J/e q d) so 
as the type-I condensation condition n >n ( iJ is fulfilled. However at = 2.52n t (at 
J=4.24kA/cm 2 ) the coherence length reduces down to the size of the cavity. Thus only in the 
narrow range nt) < n < nir , the SR pulse peak power exhibits quadratic growth (|24l) . At 
high injection levels, the peak power is oc (n — n t ) 3 ^ 2 [see the asymptotic expression (1351)1. 
This feature is perfectly reproduced by the numerical model in Fig. |U 

Fig. [5] shows the analytical model predictions for the SR output power and pulse width. 
Although the SR pulse power and pulse width calculated from the numerical model in Figj4] 
exhibit large amplitude and timing instabilities, they are in agreement with the predictions of 
analytical model. Thus at the current density 30 kA/cm 2 , both models predict the SR pulse 
peak power of ~200 W. The pulse width 2r p of 450ps predicted by analytic model (Fig. \5\ 
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Figure 4. Numerical simulations: Peak power (a), FWHM pulse width (b) and delay time 
(c) of SR pulses emitted at the absorber section facet (blue data points and curves) and gain 
section facet (red data points and curves) in function of the pump current density. In (a) and (b), 
we introduce the eye guides (doted black curves) to simplify interpretation of the numerical 
simulations subjected to large stochastic variations. The cavity length is 800 /im, the absorber 
length is 160 ^m. 
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Figure 5. Analytical solution of GLE d!51 i: Peak power (a) and pulse width r c at 42% of 
the peak power (b) of the output SR pulses in type-I [Eq.d24lil and type-II [Eq |29l regimes in 
function of the initial carrier density no. The carrier densities are shown normalized on the 
transparency carrier density n t . The sample size is 640 /im , corresponding to the gain section 
lengths of the cavity used in numerical simulations in Figs. |2j[3]and|4] 



(b)) corresponds to the SR pulse FWHM 1.76r p of 400 fs ll26ll . This nearly reproduces the 
lower boundary estimates for the pulsewidth obtained from the numerical model simulations 
in Figj4](b) [see the eye guiding curve]. 

Figure [6] illustrates a relationship between the critical densities for type-I and type-II SR 
as a function of the overall cavity length L cav , assuming that absorber occupies 20% of the 
cavity (L = 0.8L cav ). For comparison, the threshold carrier densities in the CW regime are 
plotted at lowest (V a = 0) and highest (V a = — 1) absorption states of the absorber. In most of 
the cases, the type-II SR regime with condensate fraction occupying only a part of the cavity 
will be reached provided the cavity is made sufficiently long to prevent lasing. 

The type-I SR regime, in which the condensate of lower polaritons occupies the entire 
sample, can be achieved in short cavities. Fig|7] shows predictions of analytic model for the 
SR pulse peak power in a cavity of only 50 /im long. The first and second critical densities 
for type-I and type-II SR regimes in such short cavity are well separated from each other 
{jic) = lln t and ni^ = l^n t , respectively). Note that predictions of analytic model are in 
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Figure 6. Critical carriers densities n cr and n cr2 vs cavity length in comparison with the CW 
lasing threshold carrier densities for the cavity with relative absorber length 20% in non biased 
(Ya = 0) and highly biased (V a = — 1) states. The carrier densities are shown normalized on 
the transparency carrier density nj. 




Normalized carrier density 



Figure 7. Peak power of SR pulses in 50 length cavity. Dashed lines show the analytic 
model predictions for type-I and type-II regimes, points shows the results of numerical 
simulations, the symbols "+" and "x" distinguish SR pulses emitted at the two cavity facets. 
The carrier density is normalized at the transparency carrier density n t . 



good agreement with the results of numerical simulations in the short cavities as well (Fig|7l 
points). 

Figure [8] shows the results of numerical simulations for spatiotemporal dynamics of 
the order parameter and the output SR pulse shape in the type-I and type-II regimes of 
condensation in this cavity. 

For n = 13n t (Figs. [8] (a) and (b)) the SR pulse width is 5 ps, much longer than the 
cavity roundtrip time (1.2ps). The SR pulse builds up after 15 cavity roundtrips. The order 
parameter P reaches the maximum of 2n t , in perfect agreement with the estimates of the 
analytic model of n — n£) (see Eq.(fl4l)). The SR pulse shape well agrees with hyperbolic 
secant shape assumed in the analytic model (Fig. [8](b)). The coherence length is such that the 
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Figure 8. Numerical simulations: Spatiotemporal dynamics of the order parameter P [(a) and 
(c)] and SR pulses [(b) and (d)] in type-I [(a) and (b)] and type-II [(c) and (d)] regimes. The 
samnle size is 50 «m. The initial carrier density nn is 13m T(a) and (b)l and 28m IYc) and (d)l. 
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condensate of LP occupies the entire sample of 50 /mi long. 

In the second example plotted in Figs. [8](c) and (d), the initial carrier density n = 28n t 
well exceeds the critical density nd and condensation of LP associated with emission of SR 
pulse undergos type-II transition. All stored energy is emitted in a single-shoot SR pulse after 
3 cavity roundtrips, when the order parameter P reaches a maximum of I5n t . The peak value 
of the order parameter during condensation of LPs is very close to predictions of the analytic 
model n — n& = 14n t . Small difference can be attributed to the ringing in the output SR 
pulse envelope, which is not taken into account by analytic model. 

The numerical simulations exhibit large amplitude and timing instabilities of the output 
SR pulses. The macroscopically large fluctuations are the characteristic features of a 
condensate state [|29l l28l . In this sense, the SR is a particularly interesting quantum optics 
phenomenon because the quantum fluctuations manifest themselves macroscopically, in the 
time and energy domains. Of course, these features cannot be reproduced by analytic model 
based on Ginzburg-Landau master equation (TTST) . 

The delay time statistics and timing jitter of SR pulses in atomic (molecular) ensembles 
have been thoroughly studied previously by both classical and quantum mechanical 
approaches. In some conditions, the standard deviation of the SR delay time was found to 
be 10-12 % 11301 |3TI . However, the results of previous studies can hardly be applied to our 
case of semiconductor laser structures. The reason for this is as follows. 

One of the fundamental differences between SR in atomic (or molecular) gases and SR 
in semiconductors considered here consists in the opposite relation between the decoherence 
time T 2 of the medium and the characteristic width of superradiant pulse r c . Indeed, r c <C T 2 
in the classical case of SR emission in a gas medium. The SR pulse width r c is so large 
(up to few hundreds of nanoseconds) that the coherent domain L* occupies entire sample. 
All quantum oscillators contribute coherently to the field of the SR pulse, leading to the 
superposition / oc \Ei\ 2 ) = (\Ei\ 2 ). From this point of view, the SR emission in gases 
is similar to the type-I SR regime in semiconductors (Sec. 13.21) . However, in semiconductors, 
the SR pulsewidth r c exceeds the natural decoherence time T 2 (~ 100 fs). A more gentle 
growth of the output peak power with the initial carrier density in the type-II SR regime (|40"1) 
has been ascribed to a partial coherence in the ensemble of e-h pairs during SR emission (Sec. 
13.31) . Extending these considerations, we shall conclude that some part of the initially stored 
energy in the system as non-equilibrium electron-hole pairs of density no is spent to sustain 
the coherence of SR pulse against the dephasing processes. Thus in Fig. |2l the SR pulses with 
ringing evolve survive during a remarkably long time as compared to the ultrafast intraband 
dephasing processes, while the invariant relationship (fl4l) is not perfectly fulfilled. 

The predictions of our model for the SR pulse width being longer than the inherent 
dephasing time of carriers in the InGaN/GaN QWs is in agreement with the previous 
experimental studies of SR in GaAs devices ifTTl [T71 . where the following considerations 
have been used [fT2|. MacGillivray and Feld discussed the effect of high optical gain on 
dephasing and relaxation times of the system and arrived to conclusion that in the presence 
of the high optical gain the system can eventually dephase but with characteristic time 
T 2 eff ~ T 2 g (n — n t )YL/v g . In conditions of our numerical experiment, this effective 
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dephasing time is 4 ps. The macroscopic polarization of e-h system eventually dephases, 
but the effective dephasing time of the whole system is increased. 

5. Conclusion 

In summary, we present the analytical model of SR emission and condensation of ID 
longitudinal polaritons in semiconductor laser cavities. The master equation for the LP 
condensation and SR emission is shown to be analytically similar to the Ginzburg-Landau 
equation or Gross-Pitaevskii equation. The effective mass of condensing lower polaritons is 
very low, 10 _11 m , giving a good reasoning for the thermal de Broglie wavelength criteria 
of macroscopic quantum degeneracy do not appear in our GLE for SR. We predicted two 
regimes of SR emission and polariton condensation. We show that analytic model predictions 
well agree with the results of the numerical simulations based on the traveling-wave Maxwell- 
Bloch equations. 

The effective mass of longitudinal ID polaritons in edge-emitting laser cavities is many 
order of magnitude smaller as compare to effective masses of 2D exciton-polaritons in 
semiconductor microcavity systems, rendering ID edge-emitting laser structures particularly 
attractive for experimental observations of macroscopic coherent states at room temperatures. 

The work on incorporation of the finite interband relaxation time of carriers into the 
model is ongoing and will be published elsewhere. Other intersting development is related 
to the following. The GLE equation for BCS transition in electron system contains vector 
potential of the field, so as the critical condensation density is a function of the applied 
magnetic field. Similar effect in polariton system (fl~5l) is under investigation now. 

Acknowledgments 

This research is supported by the EC Seventh Framework Programme FP7/2007-2013 under 
the Grant Agreement n 238556 (FEMTOBLUE) 

References 

[1] S. Christopoulos et al, Phys. Rev. Lett. 98, 126405 (2007). 

[2] D.L. Boiko, PIERS online 4, 831 (2008). 

[3] P.P. Vasil'ev, Rep. Prog. Phys. 72, 076501 (2009). 

[4] R. H. Dicke, Phys. Rev. 93, 99 (1954). 

[5] N. Skribanowitz et al, Phys. Rev. Lett. 30, 309 (1973). 

[6] M. F. H. Schuurmans, Q. H. F. Vrehen and D. Polder, Adv. At.Mol. Phys. 17, 167 (1981). 

[7] A. V. Andreev, V. I. Yemel'yanov, and Yu. A. Il'inskii, Cooperative Effects in Optics: Superradiance and 
Phase Transitions (Bristol: Institute of Physics Publishing, 1993). 

[8] M. Gross and S. Haroche, Phys. Rep. 93, 301 (1982). 

[9] L. I. Men'shikov, Phys. Usp. 42, 107 (1999). 

[10] P.P. Vasil'ev et al, Phys. Rev. B 64, 195209 (2001). 

[11] P. P. Vasil'ev, Phys. Status Solidi B 241, 1251 (2004). 

[12] J.C. MacGillivray, and M.S. Feld, Phys Rev A 14, 1 169 (1976). 



Longitudinal polariton condensation and superradiant emission 



23 



[13] B.M. Pimentel and A.H. Zimerman, "Phase transition in a generalized Dicke model of superradiance," 

Phys. Lett. 53A, 200(1975) 
[14] K. Hepp and E. H. Lieb, "On the superradiant phase transition for molecules in a quantized radiation field: 

the Dicke Maser model," Annals of Phys. 76, 360(1973). 
[15] Y. K. Wang and F. T. Hioe, "Phase Transition in the Dicke Model of Snperradiance," Phys. Rev. A 7, 831 

(1973) 

[16] R. Graham and H. Haken, "Laserlight-first example of a second-order phase transition far away from 

thermal equilibrium," Z. Physik 237, 31 (1970) 
[17] P. P. Vasil'ev, Quantum Electron. 29, 842 (1999). 
[18] W. G. Scheibenzuber et al, J. Appl. Phys. 109, 093106 (201 1). 
[19] J. Dorsaz et al, Appl. Phys. Lett. 98, 191 1 15 (201 1). 

[20] IV. Smetanin, P.P. Vasil'ev, and D.L. Boiko, Opt. Express 19, 171 14 (201 1). 
[21] M. Kuramoto et al, Appl. Phys. Lett. 96, 051 102 (2010). 

[22] E.T. Jaynes, F.W. Cummings, "Comparison of quantum and semiclassical radiation theories with 

application to the beam maser". Proc. IEEE 51, 89109 (1963) 
[23] M. Tavis, and F.W. Cummings, "Exact Solution for an N-Molecule-Radiation-Field Hamiltonian," Phys. 

Rev. 170, 379(1968). 

[24] F. T. Arecchi and E. Courtens, "Cooperative phenomena in resonant electromagnetic propagation," Phys. 
Rev. A 2, 1730(1970) 

[25] F. T. Arecchi and R. Bonifacio, "Theory of optical maser amplifiers," IEEE Joourn. Quantum Electron. 
QE-1, 169(1965) 

[26] In definition of the coherence time of the pulse, there is an element of ambiguity related to the intensity 
level at which the pulse width is measured. Thus 2r p is the with of the pulse A\(t— z/v g ) at the level of 
l/cosh 2 (l)«42%, while the FWHM is 2 \n.(\+y/2)T p w 1.76t p . The difference between two values is 
not fundamental. 

[27] Practical realization of SR emission from semiconductor laser cavity requires an intracavity absorber. In 
that case, the size of the sample treated by the theory is defined by the length of the gain section. 

[28] B. L. Altshuler, P. A. Lee and R. A. Webb, Mesoscopic Phenomena in Solids, ed V. M. Agranovich and A. 
A. Maradudin (Amsterdam, 1991) 

[29] L. V. Butov et al, Phys. Rev. Lett. 73, 304 (1994). 

[30] F. Haake et al, Phys. Rev. A 23, 1322 (1981). 

[31] E. A. Watson et al, Phys. Rev A 27, 1427(1983). 



